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1. Introduction 

This research is concerned with adaptive learning ' control systems for 
aircraft* The research to date, has led to the development of a learning 
control system which blends the gain scheduling and adaptive control into a 
single learning system that has the advantages of both* An important feature 
of the developed learning control system is its capability to adjust the gain 
schedule in a prescribed manner to account for changing aircraft operating 
characteristics* Furthermore, if tests performed by the criteria of the learn- 
ing system preclude any possible change in the gain schedule, then the overall 
system becomes an ordinary gain scheduling system. 

The research accomplished to date is presented in the following sections 
First, a statement of the problem is given and then the development of the 
learning control system is described. Two examples are also discussed. The 
results of the first example were presented at Langley on July 31, 1975* 

2* Statement of Problem 

Let the aircraft motion be described by 

£ » A (v) . x + B (y) • r (1) 

“P P *• -p p — 

where x is the state vector and r is the input vector. The aircraft is con- 

-p “ 

trolled using a feedforward gain matrix G and a feedback gain matrix K. The 
augmented system is therefore characterized by 

2p - [yx> + y*> • k ] ♦ [yx>- G ] i <2> 

The set of linear equations in Eg. (1) is obtained by linearizing the 

nonlinear equations of notion in the neighborhood of the operating point v. 

After this is done one obtains data points for the elements of the matrices 

A (^) and B (^) for descrete values of Let ^ £ Y, where Y is the vector 

P P 
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subspace of admissible operating conditions. Since we are interested in the 
behavior of the aircraft through the continuum of the subspace Y, we inter- 
polate the data points and therefore the elements of the matrices A^(^) and 
Bp(^) become functions of the vector Our next step is to obtain a better 
approximation of the parameters than this initial interpoation has given and 
also be able to detect changes in the elements of matrices A^(^) and B p (^) 
so as to control the aircraft effectively by updating the matrices K and G. 

For this purpose we developed an adaptive learning control system or briefly 
a learning control system which is described in the next sections. 

3. The Learning Control System 

A block diagram illustrating the functional organization of the learn- 
ing control system is depicted in Fig. 1. One of the features of this system 
is its capability to adjust the feedforward and feedback gains in a prescribed 
and learned manner to account for changing aircraft operating characteristics. 
As shown in Fig. l. t the learning system consists o£, three basic subsystems: 

1. The information acquisition subsystem, 2. The learning algorithm sub- 
system and 3. The memory and control process subsystem. The tasks of each 
subsystem and their mathematical development are described in the next three 
sections • 


3.1. The Information Acquisition Subsystem 

The information acquisition subsystem identifies the values of the 

elements of A (y) and B (y) matrices. Several techniques have been formulated — 
P *" P *■ 

all were based on the second method of Liapunov to insure convergence of the 
identification process. The technique that produced the best results is de- 
scribed in this section. In this technique the plant was represented by a model 
of the form 


x 
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F)x 
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( 3 ) 


where F is a stable matrix and the model has the same dimensionality as the 
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plant. The adaptation error is defined as 


e ■ x 




( 4 ) 


and the error differential equation describing the adaptation error is obtained 

as 
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where 
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Vectors and d^ are constant for all i, and u^, w^ are vectors whose compo- 
nents are the misallignments of the parameters of the i-th row. 

An appropriate Liapunov function for Eq. (5) is 
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where H, and are symmetric positive definite matrices with constant 

elements. Applying Liapunov's stability criterion to V and its time deriva- 

tive V, one derives a set of controller equations which when related to the 

model matrices A and B yields 
m m 
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4. 

The performance of the information acquisition subsystem is monitored by 
a convergence criterion. This criterion is defined presently by 

m V(t)ACO) £ (11) 

where ^ is prescribed by the designer according to his desired accuracy 
of identification for a prescribed length of time while the system is sub- 
jected to sufficient excitation. Other performance criteria are also being 
studied. 

# 

3.2. The Learning Algorithm Subsystem 

The identified instantaneous values of the parameters affecting the 
motion of the aircraft are fitted to predetermined analytical expressions 
describing the behavior of the parameters over the subspace Y. The predeter- 
mined analytical expressions were obtained by interpolating the apriori avail- 
able data and representing it by polynomials. In general, any set of linearly 
independent functions can be used for this purpose as long as they span the 
space of the parameter functions. 

The elements of the model matrices A and B are formed into an n*l 

m m 

dimensional vector Pj(^) given by 

- H(^)-c + v (12) 

Vector c is to be updated sequentially upon receiving new information (P , y ) . 

* 

This type of learning falls into the category of stochastic approximation 
method. There are a few ways of updating c; a survey is given in reference 
( 1 ) . Examining the different algorithms for updating vector c one observes 
that there is a tradeoff between computation complexity and rate of conver- 
gence. In our case, we used the least square error algorithm which converges 
rapidly and is of the form 
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where 


is the covariance matrix of the error vector 


I,lc 


-Vic 


(15) 


For Cq we use the coefficients of the interpolation polynomials over the 
apriori available data and Pq can be chosen any positive definite matrix. 

The choice of Pq influences considerably the rate of convergence of the 
algorithm. 

A necessary test for the learning algorithm must be to continually 
evaluate the validity of the information in vector C^. This test is per- 
formed by the confidence criterion 

0 k+ l Oc+1) 6 (16) 


where 0 fc+1 « 

The quantity 
square error 
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6 is fixed by the designer, and is the maximum tolerable mean 


3.3. The Memory and Control Process Subsystem 

After passing the confidence criterion, vector C is used to compute 

the elements of matrices A and B . The stored values of this vector (in a 

m m 

dynamic menory) are used to compute the elements of matrices A^(y) and y * 
The gain matrices K, G are then computed using the equations (6) and (7) by 
letting u^ » 0, w^ o 0. We therefore have 


6 




(17) 
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(18) 


After the gains are computed , vector C updates the previous value in the 
dynamic memory. 


4. Examples 

Two examples were studied in order to verify the validity of the learn- 
ing control system and be able to check out the programs used in the process. 

The first example illustrates the learning of three parameter curves 
for a second order representation of the longitudinal dynamics of an aircraft. 
The state variables for this system are pitch rate q and pitch attitude 0. The 
independent variable for the parameter curves is 0. The three parameter curves 

i 

to be learned are M^(0) , D^Ce) and 0^(0) — (moment, damping and control ef- 
fectiveness) . The results were presented at NASA, Langley Research Center on 
July 31, 1975. 


The second example illustrates the learning of 10 parameter curves for 
a fourth order representation of the longitudinal dynamics of the F-8 aircraft. 
The state variables of the system are pitch rate q, variational speed V, angle 
of attack , and pitch attitude 0. The parameter curves are functions of the 
mach number M and of the altitude. The equations of motion are linearized at 
selected flight conditions and can be written as 
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where 
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and 6 is a linear combination of the states and the pilot's input. The 

a^'s and b^'s are given in reference (2). In this example we used data 

from reference (2) and considered the four wing-down (CO) configurations at 

sea level; hence the 10 unknown elements of matrices A, and B, were 

long long 

functions of the mach number M. We did not use the data point at M - 1 since 
it was discontinuous with the rest of the data and it would amount for non 
valid predetermined functional representations. 

Our first step was to interpolate polynomials through the data in the 
least square sense and determine a functional representation of the system 
parameters with respect to M. After this was done we had a vector of coeffi- 
cients (C ) for the model. Then we chose a different set of coefficients 
in 

(Cp) according to which we computed the plant's unknown parameters for a given 
mach number. We simulated a flight at sea level at three different mach num- 
bers and the learning system reproduced the parameter curves (a^ & b curves) , 
An interesting observation is that the parameter curves that multiply 
the state variable V are learned in an easier fashion than the other parameters 


! 


B. 


since the response of V to a given input is sensitive to these parameters. 

To overcome the difficulty that the state variables q and of are not sensitive 
to their corresponding parameters we had to give a number of inputs at each 
operating condition so as to achieve approximate identification and therefore 
good learning. Note that the a^'s and b^*s discussed in this example are 
related to the stability derivatives for the P-8 aircraft as shown on pg. 3 
in reference (2). 

5. References 

1. Gura, I. A., "An Algebraic Solution of the State Estimation Problem," 

. AIAA Journal, Vol. 7, No. 7, July 1969, pp. 1242-1247. 

2. Gera, J., "Linear Equations of Motion for F-8 DFBW Airplane at Selected 
Flight Conditions," NASA, Langley Research Center, F-8 Digital Fly-By - 
ftire Internal Document, Report No. 010-74. 










